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Abstract 

Lamb has identified a certain class of moving space curves with soliton 
equations. We show that there are two other classes of curve evolution that 
may be so identified. Hence three distinct classes of curve evolution are asso- 
ciated with a given integrable equation. The nonlinear Schrodinger equation 
is used to illustrate this explicitly. 
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1 Introduction 



Over two decades ago, Lamb P| presented a formalism which showed that 
certain special types of motion of space curves can be mapped to com- 
pletely integrable, soliton-supporting |J, nonlinear partial differential equa- 
tions (NLPDE) such as the nonlinear Schrodinger (NLS) equation, the sine- 
Gordon equation, the Hirota equation, etc., indicating that the corresponding 
curve motions are also integrable. This formalism arose as an extension of 
Hasimoto's earlier work ||, which had established a connection between the 
equation of motion of a vortex filament regarded as a moving space curve, 
and the NLS. Since all soliton equations |2| possess similar characteristics 
such a Lax pair, solvability by the inverse scattering transform method and 
an infinite number of constants of motion, the above mapping to the NLS 
provided a motivation for a general formulation to investigate the interesting 
curve evolutions that get associated with soliton-supporting equations. 

Central to the Lamb formulation is the introduction of the following com- 
plex function if), called the Hasimoto function, defined || as 

if>(s,u) = Kexp[i J rds). (1.1) 

Here s is the arc length of the curve and u the time. K(s,u) and r(s, u) 
denote, respectively, the curvature and torsion of the moving curve. It is this 
function if) that satisfies various integrable equations in Lamb's work, and is 
used frequently in the study of various aspects of curve-dynamics j|. 

It is clear that by comparing the functional form of if) (Eq. (1.1) with 
a soliton solution of a given integrable equation for if), one can identify 
K and r, and hence obtain the associated moving curve parameters that 
correspond to an integrable, shape-preserving curve motion. This in turn 
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unravels a certain special geometric structure of the given soliton-bearing 
NLPDE. In this Letter, we show that two other complex functions defined 
by $(s,w) = rexp [i J Kds] and x = {K — it), also arise from the basic 
curve evolution equations, just as naturally as the Hasimoto function tp does. 
We demonstrate that these can also satisfy various soliton equations. This in 
turn leads to the result that each such integrable equation is associated with 
not just one, but three distinct classes of space curve motion, and therefore 
has a much richer geometric structure than hitherto envisaged. We illustrate 
this by using the NLS as an example. 

2 Moving space curves 

A space curve embedded in three-dimensions may be described using the 
usual Frenet-Serret equations || for the orthonormal triad of unit vectors 
made up of the tangent t, normal n and the binormal b : 

t s = Kn ; n s = -Kt + rb ; b s = -rn. (2.1) 

Here, s stands for the arc length of the curve. The subscript s denotes (d/ds). 
K and r are the curvature and torsion of the curve. 

If the curve moves with time it, then all quantities in Eq. (2.1) become 
functions of both s and u. To describe the time evolution of the triad (t, n, b) 
we write down the following set of equations ||: 

t u = gn + hb; n„ = -gt + r b; b u = -ht - r n. (2.2) 

Here, the subscript u denotes (d/du). As is clear, the parametes g, h and r 
which determine the motion of the curve are also functions of both s and u. 
On requiring the compatibility conditions 

t us , n sn n us , b su b us , (2.3) 
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a short calculation using Eqs. (2.1) and (2.2) leads to 



K u = (9s ~rh) ■ r u = (T ) s + Kh; h s = (Kr -rg). (2.4) 

3 Formulation (I) using the Hasimoto func- 
tion tp. 

In Lamb's formulation [JEf, which will be referred to hereafter as formulation 
(I), the second and third equations of the set (2.1) are combined to yield 

(n + ib) s + ir(n + ib) = -Kt. (3.1) 

This immediately suggests the definition of a certain complex vector 

N = (n + zb)exp[i [rds]. (3.2) 

Differentiating Eq. (3.2) with respect to s and using Eq. (3.1), we get 

N s = -Kexp[i J Tds]t. (3.3) 

Thus the Hasimoto function ip (Eq. (1.1)) appears in a natural fashion in 
the above equation. Using the definitions of N and ip given in Eqs. (3.2) and 
(1.1) respectively, Eqs. (2.1) can be written in the form 

t s = ^*N + ^N*) ; N s = (3.4) 

Next, Eqs. (2.2) take on the form 
1 



where 
7i = 



2 ( 7l *N + 7l N*) ; N u = zi^N + 7 it, (3.5) 



-(g + ih)exp[i J rds] ; R\ = (J r u ds — r ) = J Khds. (3.6) 



Here we have used the second equation in the set (2.4) to write the last 
equality in (3.6). We use the subscript 1 here, to denote formulation (I). 
From Eqs. (3.4) and (3.5), setting N su = N us , and equating the coefficients 
of t and N, respectively, we get 

ipu + lis - iRiip = 0; (3.7) 

Ru = ^(Ti^-TiVO- (3-8) 

Interestingly, as noted by Lamb 0, the structure of Eqs. (3.7) and (3.8) 
which arose from compatibility conditions on curve evolution, suggests a 
possible relationship with soliton-bearing equations, via the Ablowitz-Kaup- 
Newell-Segur (AKNS) formalism . This is seen as follows : It is well 
known that for a class of soliton-bearing equations with a function q(s, u) as 
the dependent variable, the Lax pair L and M in the AKNS formalism are 
given by : 

^ = ( B t U,C l)v = M y (3.10) 

8u \-B*(s,u,() -A(s,u,()) y y V ; 

Here, the eigenfunction y is the column vector (yi y<i) T and £ is the eigenvalue 
parameter. Requiring y su = y us , Eqs. (3.9) and (3.10) lead to the following 
AKNS compatibility conditions 0: 

q u = 2Aq + B s + 2iCB (3.11) 

A s = (Bq*-B*q) (3.12) 

Equations (3.11) and (3.12) are identical in form to Eqs. (3.7) and (3.8) 
provided the following identifications are made : 

q = ?P/2; A = iR 1 - B = - Tl /2; Q = 0. (3.13) 
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Thus the curve evolution equations (2.1) and (2.2) imply AKNS equations, 
with C = 0. Now, 71 and R\ are given in Eqs. (3.6). It was shown by Lamb 
[[!]] that for appropriate choices of 71 as a function of ip and its derivatives, Ri 
can be found from Eq. (3.8), and substituted in Eq. (3.7) to yield integrable 
equations for ip. 

Alternatively, once the quantities 71 and R\ are identified so that Eqs. 
(3.7) and (3.8) take on the form of an integrable NLPDE for ip, it is possible 
to find the corresponding Lax pair, by working directly with equations (3.4) 
for t s and N s , along with equations (3.5) for t u and N u : This is done [|l| by 
considering the constraint tf + nf + bf = 1, where the subscript I = 1, 2, 3 is 
used to represent the three components of the vectors concerned. One then 
defines a complex function 

= (n, + ik)/(l - t{) = (1 + ti)/{ni - ik), (3.14) 
(2) (3) 

along with its two other counterparts /j and , obtained by cyclically 
changing t, n and b in Eq. (3.14). It can then be shown that all the three func- 
tions fi, a = 1, 2, 3 satisfy appropriate Riccati equations. The correspond- 
ing Lax pair in each case can then be obtained by setting = y^/yi, 
and identifying the corresponding eigenfunction to be the column vector 
{V\^ 2/2^) 7 '■ ^ short calculation shows that only one of the Lax pairs thus 
obtained has the AKNS form of Eqs. (3.9) and (3.10), with entries as in 
(3.13). 

We paranthetically remark that it is possible to introduce a non-zero 
eigenvalue parameter ( in this formalism, by means of a suitable gauge trans- 
formation of ij} in Eq. (3.7). This in turn can be related to a certain gauge 
freedom in the choice of the orientation of the axes in the plane per- 
pendicular to the tangent of the space curve. The details of this and its 
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ramifications will be reported when completed. 



4 New formulations (II) and (III) using func- 
tions <E> and % 

As already noted, in Lamb's formulation or formulation (I), one proceeds by 
first combining the second and third equations of the Frenet-Serret set of 
equations (2.1), leading to the appearance of ip. In this section, we consider 
the other two possibilities: Formulation (II), that combines the first and 
second equations of (2.1), and formulation (III) that combines the first and 
third equations. As we shall see, these formulations lead to the appearance 
of two other complex functions $ and x respectively, in a natural fashion. 
Formulation (II): Combining the first two equations in Eqs. (2.1), we get 

(n-it) s + iK(n-it) = rb (4.1) 

This suggests the definition of a second complex vector 

M = (n-it) exp[i J Kds] (4.2) 

Differentiating Eq. (4.2) with respect to s and using Eq. (4.1), we get 
M s = r exp[i / Kds]b. Thus a second complex function 

$(s,w) = rexpfi J Kds] (4.3) 

appears in a natural fashion in this case. Using the definitions of M and 
$ given in Eqs. (4.2) and (4.3) in the basic equations (2.1) and (2.2), and 
repeating the steps of Sec. 3 , we obtain the following counterparts of Eqs. 
(3.4), (3.5) and (3.6): 

M s = $b ; b s = -i($*M + $M*), (4.4) 
7 



M tt = ii2 2 M- 7 2b ; K = i( 72 *M + 72 M*) (4.5) 

where 

72 = — (to — ih)exp[i J Kds] ; R 2 — (J K u ds — g) = — J rhds. (4.6) 

We have used the first equation in the set (2.4) to write the last equality in 
(4.6). The subscript 2 is used on 7 and R to indicate that these correspond 
to formulation (II). From Eqs.(4.4) and (4.5), setting M su = M us , and 
equating the coefficients of b and M, respectively, we get 

$ u + 72,-^2$ = 0; (4.7) 

R2s = |(72$* -72*<*>)- (4-8) 

The structure of Eqs. (4.7) and (4.8) is the same as that obtained in 
Lamb's formulation (see Eqs. (3.7) and (3.8)). Thus using the same steps as 
in the last section, these can also be cast in the form of AKNS compatibility 
conditions (3.11) and (3.12), with the following new identifications: 

g = $/2; A = tR 2 ; B = - l2 /2- ( = 0. (4.9) 

This shows that, just as one does in Lamb's formulation, here also we could 
take choices of 72, as an appropriate function of $ and its derivatives, find 
R 2 from (4.8), and substitute these expressions in Eq. (4.7) to obtain some 
of the well known integrable equations. 

Formulation (III): Combining the first and third equations of (2.1), we get 

(t + ib) s = (K-ir)n (4.10) 
This suggests the definition of a third complex vector 

p = (t + ib) (4.11) 
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This leads to P s = (K — ir)n. Thus a third complex function 

X (s,u) = (K-ir) (4.12) 

appears in this case. Using the above definitions of P and x given above in 
the basic equations (2.1) and (2.2), and proceeding the same way as in the 
previous formulations, we get 

P s = X n ; n s = -±( X *P + *P*), (4.13) 

P u = tR 3 P- l3 n ; n u = I( 7 *p + 73 p*) (4.14) 

where 

7s =-(g-iT ) ; R 3 = -h = - J(Kr -rg)ds. (4.15) 

Here we have used the last equation in the set (2.4) to write the last equality 
above, and the subscript 3 corresponds to formulation (III). /.From Eqs.(4.13) 
and (4.14), setting P SM = P us , and equating the coefficients of n and P, 
respectively, we get 

Xu + lss-iRs X = 0; (4.16) 
Rss = '-(is X*-7 3 * X)- (4-17) 

Once again, Eqs. (4.16) and (4.17) have the same form as Lamb's result 
(Eqs. (3.7) and (3.8)), as well as the AKNS compatibility conditions (3.11) 
and (3.12). Here, the identifications are 

q = x /2; A = iR 3 ; B = - 73 /2; C = 0. (4.18) 

Thus using the same reasoning as before, for suitable choices of 73 as functions 
of x an d its derivatives, Eq. (4.16) for x can take the form of known integrable 
equations. 
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Now, as is clearly seen from Eqs. (3.6), (4.6) and (4.15), the parame- 
ters j n , n — 1,2,3 that arise in the three formulations, correspond to three 
complex functions involving different combinations of the curve evolution 
parameters g, h and r . Further, the corresponding complex functions i/j, $ 
and x that appear in the formulations are also different functions of K and 
r. (See Eqs. (1.1), (4.3) and (4.12) respectively.) Thus it is clear that the 
three formulations describe three distinct classes of curve motion. In the next 
section, we discuss an application to show this explicitly. 

5 An illustrative example: The nonlinear 
Schrodinger equation 

Next, we consider the application of our results to the NLS. This exam- 
ple is appropriate not only because the NLS was one of the first integrable 
equations to appear in the context of curve motion ||, but also because it 
has applications 0] in various fields such as vortex filament motion, optical 
solitons, magnetic chain dynamics etc. 
Formulation (I): As Lamb has shown, the choice 

7i = -tys, (5.1) 

when used in Eqs. (3.8) leads to 

Ri = \\1>\ 2 . (5-2) 
Substituting this in Eq(3.7) yields the NLS: 

iip u + ipss + l\ip\ 2 '4> = 0- (5.3) 

Using Eq. (1.1) in Eq. (5.1) and equating the resulting expression to 71 
defined in (3.6), we get the following curve evolution parameters g, h and tq 
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for this case: 

g = -Kt ■ h = K s ■ r = (h s + rg)/K = (K ss /K)-t 2 (5.4) 

On the other hand, substituting the above expressions for g and h in the 
first equation of (2.2) leads to t u = gn + hh = —Km + K s h. By using 
Eqs. (2.1) and (2.2), it can be easily verified that the above equation implies 
t u — t x t ss . This is identical in form to the Landau-Lifshitz (LL) equation 
S u = S x S ss for the time evolution of the spin vector S in a continuous 
Heisenberg ferromagnetic spin chain. In this formulation, t, the tangent to 
the space curve gets identified with S. The connection between the NLS and 
the LL equation is well known [Qj by now. 

Formulation (II): Here, the choice 72 = — i& s when used in Eqs. (4.8) leads 
to R 2 = ||$| 2 - Using this in Eq. (4.7) yields the following NLS for $: 

2$ u + $ ss + i|$| 2 $ = 0, (5.5) 

which is identical in form to Eq. (5.3) obtained in formulation (I). However, 
$ = r exp[i / Kds] here. Using this to find the above expression for 72, 
and equating it to the expression for 72 defined in Eq. (4.6), we obtain the 
following curve evolution parameters: 

h = -r a ; r = -K r ; g = {K r -h s )/r = (r„/r)- K 2 . (5.6) 

Thus the above parameters g, h and r which describe the curve motion in this 
formulation are quite different from those of formulation (I) (See Eqs. (5.4)). 
In addition, using the above expressions for h and To in the third equation 
of Eqs. (2.2), a short calculation shows that for this curve, b u = — b x b ss . 
Thus, in contrast to formulation (I), it is the binomial of the moving curve 
that satisfies the LL equation in this case. 
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Formulation (III): Proceeding as in the previous two cases, here the choice 
is 73 = ~ iXsi implying R 3 = ^\x\ 2 - This leads to 

iXu + Xss + ^\x\ 2 X = 0, (5.7) 

where x — {K — ir). (See Eq. (4.12)). The curve evolution in this case is 
described by the three parameters 

9 = r s ; h = {1/2){K 2 + r 2 ) ; r = - K s . (5.8) 

Thus by direct comparison of Eqs. (5.8), (5.4) and (5.6), we see that 
the above expressions for g, h and r , which describe the geometry of curve 
evolution in this case, are quite different from those obtained in formulations 
(I) and (II). Further, using the parameters given in Eq. (5.8), it is possible to 
show that the normal to this moving curve satisfies the LL equation in this 
case, i.e., n u = — n x n ss . We paranthetically remark that the correspondence 
of this equation to the NLS equation has been noted by Hirota in a different 



context |]I0 |. 



To demonstrate that the three formulations lead to three distinct space 
curve evolutions, each with a different set K, r, g, h and r , we consider the 
simple example of a one-soliton solution q of the NLS equation: 



q = a sech £ exp i V e (s — V c u)/2. (5.9) 

Here V e and V c denote, respectively, the envelope velocity and carrier 
velocity of the soliton. The amplitude a = [V e {V e — 2V C )}^ and £ = (s — 
V e u)(a /2). Also, V e (V e —2V c ) > 0. Note that any two of the three parameters 
a , V e and V c can be taken to be independent. We have taken the case 
V c — for illustration, and obtained the expressions for K and r for this 
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solution, in the three formulations, by using ip, $ and x defined in Eqs. 
(1.1), (4.3) and (4.12) respectively. These have been entered in Table 1. 
Using these, the corresponding time-evolution parameters g, h and Tq for the 
three formulations are computed from Eqs. (5.4), (5.6) and (5.8) respectively. 
These are given in Table 2. 

New geometries associated with a soliton solution of NLS: We now 

present the geometrical consequences that arise from the three formulations, 
at the curve level. Let the moving curve be described by the vector r(s,u). 
Thus the tangent t appearing in the Frenet-Serret equations (2.1) is defined 
as t = r s (s,u). According to the fundamental theorem of curves ||, smooth 
functions K(> 0) and r define a curve r(s,u) uniquely, modulo orientation 
in space. Now, from the expressions for K and r given in Table 1, we see 
that formulation (I) yields a moving space curve with constant torsion r but 
a space-time varying curvature K. On the other hand, the moving curve ob- 
tained using formulation (II) has a constant K but varying r. In formulation 
(III), both K and r are space-time dependent. Since the three moving curves 
have different (K, r) parameters, our results explicitly illustrate that they are 
indeed geometrically distinct. Further, since they are all gauge equivalent to 
the NLS, they correspond to integrable curve evolution. 

Now, any moving curve arising from formulation (I) is well known to be 
the solution of the following "localized induction (LI) equation" || for a 
vortex filament: r u = r s x r ss . This is because the LI equation can be 
easily shown to imply the LL equation for the tangent t that appears in 
formulation (I). In contrast, in formulations (II) and (HI), the LL equation 
is satisfied by b and n respectively, and hence the corresponding r(s, u) does 
not obey the LI equation any more. Thus the curves arising from the two 
latter formulations are distinct from that associated with the LI equation. It 



13 



is of course, not at all obvious whether analogs of the LI equation (i.e., a PDE 
at the curve level) can be written down for these two formulations. However, 
irrespective of this, it turns out that there is a procedure to construct the 
corresponding moving curve r(s, u) itself, for each of these formulations. The 
results of this lengthy calculation will be presented elsewhere. 

Clearly, the dynamics of these new curves which get associated with the 
LL equations for n and b respectively, would also be form-preserving curves, 
just as Hasimoto's vortex filament was. These may find applications in fields 
other than fluid dynamics. For instance, a phenomenological modelling of the 
motion of an interface (regarded as a curve) could be one possible application. 

6 Conclusions 

In this Letter, we have found new connections between moving space curves 
and solitons using two natural extensions of Lamb's formulation. Our unified 
procedure demonstrates that a specific evolution of an integrable equation 
can actually be associated with three different classes of space curve motion, 
the class discussed by Lamb being one of these. As described in Sec. 3, the 
corresponding Lax pairs in the three cases can also be constructed. (See Eq. 
(3.14) and the discussion below it 0.) 

Application to the NLS is discussed. Our analysis shows that the class 
of curves that are solutions of the well known " localized induction equation" 
for a vortex filament is but one of the three classes of moving curves that can 
be associated with the NLS, i.e., the class connected to Lamb's formulation. 
There are (as we have shown in general) two more classes, and in order to 
understand these new geometric structures of the NLS more explicitly, we 
have determined the various curve parameters for each of the three moving 
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curves associated with a one-soliton solution (5.9) of the NLS and displayed 
them in Tables 1 and 2. While the curve corresponding to Lamb's class has 
a constant torsion but varying curvature, the second curve has a constant 
curvature, but varying torsion. For the third curve, neither of these param- 
eters is a constant. As this example shows, these three are clearly distinct 
curves. It would indeed be instructive to apply our results to other integrable 
equations as well, and find the three associated geometric structures. 

We conclude with a general remark on the possible connection between 
a singularity in the "velocities" (i.e., the time evolution parameters) of the 
geometric dynamics and that of the associated NLPDE. We have shown that 
the three complex functions ip = K exp [i J rds], $ = rexp [i J Kds] and 
x = (K-ir) = (K 2 + t 2 )^ exp [i J tan' 1 (r/K)ds], satisfy the same NLPDE, 
and that three different moving curves get associated with it. Therefore if a 
certain solution of the NLPDE has a singularity at a point, then it implies 
that as the curve evolves, a corresponding singularity will appear in the 
curvature K of the first moving curve, in the torsion r of the second curve, 
and in the quantity [K 2 + r 2 )^ for the third. Since these "velocities" g,h 
and To will be certain functionals of K, r and their derivatives for the three 
curves, the corresponding singularity will be reflected in these velocities as 
well, and can be studied case by case. 
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Formulation 


Solution of NLS 


K 


T 


I 


1> 


Ksech£ 


V e /2 


II 


$ 


V e /2 




III 


X 


V^sech£cosV e s 


— V e sech£sinV c s 



Table 1: Example: The curvature K and torsion r for the special soliton 
solution (Eq. (5.9)) of the nonlinear Schrodinger equation (NLS) for ip, $ 
and x i n the respective formulations (I), (II) and (III). This soliton has a 
vanishing carrier velocity and an envelope velocity V e . 



Formulation 


9 


To 


h 


I 


(-K72)sech£ 


(-K 2 /2)sech 2 e 


(-y e 72)sech^tanh^ 


II 


(-K 2 /2)sech 2 e 


(-K72)sech£ 


(K72)sech£tanh£ 


III 


(K72)sech£x 
(tanh£sinV c s — 2cosV c s) 


(K72)sechCx 
(tanh£cosV c s + 2sinV c s) 


(Vj2)sech 2 £ 



Table 2: The coresponding time evolution parameters g, Tq and h in the 
three formulations (Eqs. (5.4), (5.6) and (5.8)) for the example in Table 1. 
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